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Abstract 

We study the possibility for a semi-infinite block of linear viscoelastic ma- 
terial, in homogeneous frictional contact with a non-deformable one, to slide 
under shear via a periodic set of "self-healing pulses", i.e. a set of drifting 
slip regions separated by stick ones. We show that, contrary to existing ex- 
perimental indications, such a mode of frictional sliding is impossible for an 
interface obeying a simple local Coulomb law of solid friction. We then dis- 
cuss possible physical improvements of the friction model which might open 
the possibility of such dynamics, among which slip weakening of the friction 
coefficient, and stress the interest of developing systematic experimental in- 
vestigations of this question. 
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I. INTRODUCTION 



A few recent qualitative observations [|TJ, O on the frictional motion of sheared gels 
sliding along smooth glass surfaces point towards the existence of inhomogeneous modes of 
frictional sliding. Namely, in some limited range of values of -small- shearing rates, sliding 
seems to occur via propagation of a quasi-periodic pattern of sliding zones of finite width, 
separated by non moving regions, where the interface sticks. These "slip pulses" drift at 
velocities c ~ mm/sec, while the remote average (pulling) velocity lies in the 1 — 10fim/sec 
range. Their width is typically tens of micrometers. Analogous observations have been made 
by Brune et al on a sliding rubber foam, and by Mouwakeh et al. f§] on the elastomer 
polyurethane. 

The topology of such sliding modes is reminiscent of that of Schallamach waves ||, which 
have been documented in the case of some very compliant transparent rubbers sliding on 
smooth glass. They consist of quasi periodic zones, of width typically I ~ 100/im, with space 
periods roughly ~ 10/ j7|, where the rubber buckles, so that the two surfaces get separated 
by a distance comparable with /. These separation waves have drift velocities ~ mm/sec, 
for remote velocities ~ /zm/sec. 

However, the slip pulses in gels do not seem to be associated with any interface separation. 
In this respect, they are more comparable with the so-called "self-healing slip pulses", on 
which the attention of mechanicians has been focussing recently || , following the suggestion 
by Heaton |§ that some major seismic events may have occurred, not by quasi simultaneous 
sliding of the whole rupture zone, but via fast propagation of localized sliding zones of small 
extent. 

These observations all point towards a common question about the nature of frictional 
sliding, which can be schematized as follows. Consider two very thick blocks of solid materials 
with dissimilar elastic properties, in frictional contact along a planar interface of infinite 
lateral extent (Figure 1). This system bears a remote homogeneous compressive stress r| 2 , 
normal to the interface. Assume that, under the remote shear stress rj" 2 , the upper block (I) 
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slides towards x\ > at a remote point velocity Vo with respect to the lower one (II). Such 
motion can of course occur in a homogeneous mode, where stresses are uniform. Along the 
(homogeneous) interface, the friction law, which we assume to obey the Amontons-Coulomb 
proportionality between shear and normal loads, imposes that 

T 12 = -fd(v )T2 2 (1) 

If such is the case, given r| 2 and vo, the remote shear sliding stress is fixed. In the Coulomb 
approximation, where fine variations of the dynamic friction coefficient are ignored, fd re- 
duces to a constant. 

The question then arises of whether or not this homogeneous sliding mode is stable with 
respect to small non homogeneous perturbations of the stress and strain fields localized in the 
surface region. In other words, do deformation waves exist along a sliding frictional interface? 
If so, are they damped, or amplified, or neutral? This question has been studied extensively, 
for dissimilar linear elastic materials, with Coulomb friction, by several authors, in particular 
Weertmann |TIJ, Adams [IT], and Martins et al [I2|, whose results are synthetized in a recent 



article by Ranjith and Rice ||. They find that, when fi^ 7^ and when such interface waves 
exist [T3|], the corresponding sliding velocity field along the interface has, for a mode of 



wavelength k, the form : 

v (xi) = Vk exp [ik (x\ — ct) + a \ k | t] (2) 

Given the elastic moduli, the drift velocity c and amplification coefficient a are real positive 
constants. That is, waves drifting along (resp. against) the direction of v are amplified 
(resp. damped). Homogeneous sliding is thus linearly unstable against perturbations of all 
wavelengths. 

These results are derived under the assumption that the interface is sliding everywhere. 
As amplification proceeds, the sliding velocity necessarily vanishes at some points. This sug- 
gests that sliding might occur via a periodic set of self healing slip pulses, separated by stick 
regions. A family of such pulses has been built by Adams |]TTJ for dissimilar elastic solids. 



Their drift velocity c, which depends on the values of elastic moduli, is, roughly speaking, 
on the order of a sound velocity. So, their dynamics is controlled by inertia. However, such 
self-sustaining (stationary) dynamical patterns are singular in the following sense. Since 
perturbations of all wavelengths, however small, are amplified, any initially localized per- 
turbation gives rise to diverging oscillations at arbitrary small time : Adams's pulses have 
zero measure attractors. This so-called "ill-posedness" most likely signals that the Coulomb 
friction law misses some of the physical processes which control the fast dynamics of frac- 
ture at frictional interfaces between elastic materials, i.e. their high frequency response - a 
problem which is currently under study || . 

Slip pulses in gels or rubbers, on the contrary, are slow dynamical objects whose velocities, 
comparable with those of Schallamach waves, are much too low for inertia to be relevant. 
Their dynamics is certainly controlled by the dissipation associated with the viscoelasticity 
of these materials. 

We therefore concentrate, in this article, on the following question. Let block (I) be an 
incompressible linear viscoelastic material with, for simplicity, a single viscous relaxation 
time. It slides slowly on a smooth non-deformable material, and interface friction obeys 
a simple local Coulomb law. Under such conditions, are non inertial periodic slip pulses, 
stationary in the drifting frame, a possible mode of motion? 

In Section II we formulate the corresponding mathematical problem, and derive the form 
of its analytical solutions. We show in Section III that none of these is compatible with the 
stick conditions to be satisfied in the non moving parts of the interface. Hence, in this as 
well as in the inertial regime, a Coulomb law with a constant dynamical friction coefficient 
is incompatible with the existence of such modes of motion. We discuss in Section IV 
possible physical tracks towards improvements of the simple Coulomb model, which might 
be relevant to the problem of inhomogeneous sliding, and stress the interest of corresponding 
experimental investigations. 
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II. GENERAL FORMULATION 



We follow closely the approach of Adams Jl "I J and of Comninou and Dundurs jnj re- 
stricted to the case where (see Figure 1) block (II) (x 2 < 0) is non deformable. Block (I) is 
submitted to the uniform remote stresses r 22 < and r* 2 . It is infinitely extended along x\, 
and made of an incompressible material, with a linear viscoelastic shear response described 
by: 

t 

ri2 (t) = J dt' fl{t-t')u 12 {t') (3) 

— oo 

T12 and u\2 are the shear stress and deformation, confined to the (a; 1,2:2) plane, and we 
model the time dependent shear modulus as a single time Kelvin one, namely : 

A* (t) = H00 + (]J>o - yUoo) e~ t/T (4) 

We will moreover assume that the relaxed modulus /ioo is much smaller than the short time 
one, /i . To fix ideas, for compliant rubbers, values of /x^/ /x ~ 10~ 3 are typical. 

We want to study dynamic patterns, where (I) slides towards X\ > with the uniform 
remote velocity vq, with space period A = 2ir/k, and drift velocity c in the frame of (II). 
The corresponding form of the displacements Ui, u 2 reads : 

ui =v t+Y, D ml (x 2 ) (5) 

m>l 

u 2 = D 02 + ^ D m2 (x 2 ) e imk ^-^ (6) 

m>l 

Solving the wave propagation equation together with the condition of non separation at 
the interface : u 2 (xi, x 2 — 0, t) = 0, one obtains straightforwardly (see Appendix A), in 
the incompressible limit (Poisson coefficient v = 1/2) for the interfacial sliding velocity 
v s (rj) = dui/dt \ X2=0 : 

v s (77) = v + cRe £ B m e irm > (7) 

m>l 
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where rj = k (xi — ct) while the interfacial shear and vertical stresses read : 



712 (v) = T 12 + Re ~ iB mPm (1 + <T m ) & 

m>\ 



imr\ 



(8) 



T22 (rj) = r 2 * 2 + Re B mPm (1 - <r m ) e 



imrj 



(9) 



m>l 



with 



per 

pm 



1/2 



-Re<j m > 



(10) 



We restrict our attention to the very slow modes observed experimentally, for which, what- 
ever the frequency (mck), pc 2 << p m , so that 

,,2 



cr m — 1 - 



2p r< 



(11) 



p is the mass density of material (I), p m its (complex) elastic modulus at frequency (mck) 



Pm = P (mck) 



-iu / dt p (t) e 
Jo 



hot 



-mck 



(12) 



The unknown coefficients B m must be determined from the second boundary condition 
along the interface, which we want to describe a set of slipping regions of length 21 separated 
by sticking ones. 

We assume friction to be described by a simple local Coulomb law, with a constant 
dynamic friction coefficient / equal to the static one, that is : 
(i) Slip regions : —a + 2pn < r] < a + 2pn 



r s (v) = T12 (rj) + fr 22 (rj) = Q v s > 



(ii) Stick regions : —a + 2-p-a < i] < —a + 2 (p + 1) n 



v s = fr 22 < 7i2 < ~fr 2 2 



(13) 



(14) 



This description of interface friction calls for a few comments. Indeed, it assumes tacitly 



- as is common in contact mechanics |15] - that one can legitimately define a local and 



space-independent friction coefficient. Since solid friction results from the average effect 
of disipative flips of bistable pinned elastic units, this can be true only on a scale much 
larger than (i) the size b of the basic unit, and (ii) the scale L of interface inhomogeneities. 
The detailed analysis |TjJ of the Rice-Ruina phenomenological law of dynamic friction fl? 



has shown that b is of nanometric order. So, our assumption is justified for interfaces with 
homogeneous intimate contact. Such is indeed the case for the gels or very compliant rubbers 
which we have in mind here, as long as elastic deformations vary on scales much larger than 
nanometers - which sets a lower limit on the size of Dugdale-Barenblatt-like fracture head 
regions. 

Note, however, that the situation is different when dealing with multicontact Greenwood- 
like interfaces ||18|| . These prevail with stiff materials, such as metals, glasses or rocks, which 
are not polished down to nanometric roughness. Then, the small scale cutoff is provided 
by the average distance between contacting asperities, commonly lying in the 100/zm/sec 
range. This, in our opinion, should be kept in mind when attempting, for such interfaces, 
to regularize the above mentioned ill-posedness problem, since, on space scales < L pinning 
strength fluctuations become non negligible. 

Taking condition flTT| ) into account we set, following Comninou and Dundurs for 
the periodic function v s (f]), in —it < rj < tt : 

v s (rj) = a <\ r] |< 7r (stick) (15) 
v s (rj) = v(rj) — a < rj < a (slip) (16) 



Hence, from equation 



B m = — I div(i)e- m ^ (m>l) (17) 



TIC 



vo = ^r^v(0 (18) 

271 J —a 

Using equations (|15"D, ( |16"D together with @, @, and with the help of the relation 
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m>l 



1 00 Z 



cotg- 



x 



(19) 



(where PV designates the Cauchy principal value), one finally gets, for the interfacial "sliding 

stress" 7% 



Ts(v) = r: + ^[V(r ] )-V ] 



1 ^7 

- [ d£V(0 -fdtfi(t)- 

7T J -a J at 





COt(? 



7] — £ + c/ci 



(20) 



where we have set V(?y) = v(r))/c, and^f / dx = PV J f dx. 

Integrating by parts the last term in the r.h.s. of Eq. ([TED , the condition Eq. flTU| ) for 
frictional sliding within the slip pulses provides us with the integral equation to be satisfied 
by the interfacial reduced velocity field in (—a < r\ < a), namely : 



r; + IV(V) -V ] + ^ 

2 7T 



a 00 

j diviocotg 7 ^ + - j 



dt 



dfi 
~dt 



cotg 



r] — £ + ckt 



(21) 



with : 



I" d£V(0 = 2nV 

J —a 



(22) 



Once Eqs. (|T7D are solved for V(r}), interfacial stresses in the stick regions should be 
calculated from Eq.([L6|), Eq.(|Tl]) then providing the final condition for slip pulses to exist. 

Expression ( pip separates explicitly the instantaneous elastic shear effects (3 rd term) 
from the contribution of viscoelastic relaxation (4 th term). The second term, which derives 
from the perturbation of the normal stress r 2 2, is, for our very slow pulses, smaller than the 
integral ones by a factor (c/c s ) 2 , where c s is some sound velocity. We will therefore neglect 
it from now on. 

The cotg form of the elastic kernels results from imposing space periodicity to the pat- 
terns. Eq. (pip can be rewritten in a form more standard in fracture mechanics by setting, 
in (—71 < i] < 71") : 



V Z V(u) 

u = tq- y = tg- $(u) = — 

y 2 y y 2 K J 1 + u 2 



a kl 
a = to— = to — 
y 2 2 



(23) 



Some elementary algebra then leads, in (—a < u < a), to : 



1 + u 2 n J u — y 

with : 



Uy^L + /■*» M %»)-^ = (24) 

7 it — y J-a 1 + u 



= ttVq (25) 



and we have set : 



zJ u-y + (l + uy)tg ( 



rA7 
2 



(26) 



The singular integral equation ([24]) belongs to a class which was studied extensively 
by Mushkelishvili [fLSfl . In his terminology, the first two terms on the l.h.s. constitute 
the "dominant" part. The viscoelastic kernel k(u, y) satisfies the regularity condition : 
lim u ^ y [(u — y)k(u, y)] = 0. This entails that its plays no role in the strength of the singular- 
ities of the solutions - i.e., as is intuitively reasonable, these are ruled by the instantaneous 
elastic response of the deformable medium. 

Following reference [113, there are four families of solutions of equation (1211), each of 



which is associated with one of the basic functions, characteristic of the dominant part : 

Z et , eh (y) = (y + ar /2 (*-yy h/2 (27) 

where the indices (e^) = (+1, —1) control the convergent or divergent behavior of Z at the 
tail and head edges of the slip zone. 

One can then transform, for each family, Eq.([24]) into an equivalent non-singular Fred- 
holm integral equation. It moreover turns out that, when we specialize to the single relax- 
ation time model for n(t) (Eq.@), analytical expressions for the solutions of these equations 
can be obtained explicitly, thus allowing us to draw explicit conclusions about their existence. 

In view of the heaviness of the (otherwise straightforward) algebra involved, we will 
exemplify the method in full detail only for one of the families, namely the (H — ) one. 
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III. THE FOUR FAMILIES OF SOLUTIONS : ^-FIELDS AND EXISTENCE 

CONDITIONS 



A. The (+-) family 
The corresponding basic function 

z+-(y) 



{y + a) (28) 



\l - v) 

The implementation of Mushkelishvili's method is performed in Appendix B. For fi(t) as 
specified by Eq.(f|), the non singular equation equivalent to Eq.([24"D reads : 

A/z 2exp(2tan~ 1 u/ ckr) f a 
fi ckr(l + u 2 

with : 

2/i 1 4~ u /i ckr 7r 1 4~ it 

(30) 



*+_(«) = H + 4u) + ^ gg^g^Ml) r d z*+-{z) exp (-2tan-h/ckr) (29) 

u n ckr[\ + u 2 ) Ju v ' 



W+_ = 2 J a dz<5>(z)exp I ~ 2ta ^ Z ) ( 31 ) 



„oo p 2tan 1 */cfer l\U _ n r-a p 2tan ^/cfcr /_\T/ i n 

G{u) = -(3 / (Hf , 5T\/t— ^ -0-+P) d^T^ 77 5Tl/ 

v; M l (w-^)(l + $ 2 ) yW + a V M; i-oo (m-^)(1 + ^ 2 )V -tf-a 

(32) 

/? = (e 27r/cfeT - I)" 1 Apt — fj, — ^ (33) 
C* = 2^0/^00005^ - r^sin- (34) 

-D* = -2Vo/iooSm- - i" s *cos- (35) 

Expression (|32|) for G(m) is valid for solutions whose slip zone length 21 satisfies the condition 
21 < A/2. We assume this to hold, in accordance with experimental observations, which 
indicate values of 21 /X << 1. 
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Eq. 

(p9|) is a first order differential equation for the function /" dz<& + -(z)exp[— 2tan~ 1 z/ckr], 
which is straightforwardly solved into : 



2 An 1 

ckr fj, 1 + u 2 



$+_(» = H+-(u) + — — -^—^ I dzH + -(z)exv 



yUoo 2(tan 1 u — tan 1 z) 



(36) 



fMo ckr 

This defines a family of slip velocity fields, each of which is labelled by the four dimen- 
sionless parameters Vq = vq/c, T*//i , l/X = a/2n,cr/ X. Two of the physical parameters, vq 
and t*, are "external" : in an experiment, one imposes in general an average sliding velocity 
- hence vq is fixed - and measures r*. I, c, X, are the internal parameters of the family. This 
defines a problem of dynamical selection, namely : if sliding patterns exist, are l,c,X, and 
hence t s , uniquely defined when vq is fixed, or not? In order to clear up this important 
question, it is necessary to list the relations between them - or, alternately, the conditions 
to be satisfied by $+_ as given by Eq.(^). These are : 

(i) Two consistency conditions, expressing that the remote velocity and stresses are 
simply the k = components of the corresponding fields. This is expressed by relation 
and by an analogous equation for r s : 



<i"TT^ (37) 



IX J-oo 1 + U 

where t s (u) is related to by the first of equalities 

(ii)The interfacial stress field must also satisfy the stick inequality Eq.([L4]). One easily 
checks, with the help of Eq. (|24]) , that a divergence of $ at an edge, u = ±a, of the slip zone 
results in a diverging t s at the corresponding stick zone edge, and therefore in the violation 
of the stick condition. (Eq.(p8|)) diverges at the slip head u — a. For u — > a : 



^ / \ z +-( u ) r C*a + D* 2A/i W^ + _ ^, . . . 

with lim M ^ a 3?(n) = 0. 

A necessary condition for $ H to be acceptable is that the coefficient of in Eq.i 

vanishes, i.e., using Eqs. (|29| - |35D : 
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„ a 2 H/ + _G(a) 
r .° M 2 = -^ A/ '-dfcr- (39) 

So, for solutions of the (H — ) class, the five pattern parameters are linked by three relations. 
That is, for a given vq, this class of patterns, if they exist, form a one parameter family. We 
will comment further on this conclusion in section IV. 

Let us now come back to the "regularization condition" Eq.fl39"|). From conditions ( |T5| ) and 
(Hp, the interfacial sliding stress t s must be non positive everywhere. Hence, its w-average 



r*, must be strictly negative. 

On the other hand, the v s , and thus the $ field must, by Eq. (|i"3|) , be positive everywhere 



in the slip zone. Then definition (|3l| ) entails that W + - > 0. Finally, using Eq . (p2|) , one gets: 



G(a) = - / dV Lcfcr \ 2 ^ — < 40 



Therefore, condition fl39|) can never be satisfied. No solution of type (H — ) exists. In 
other words, viscous relaxation and pulse-pulse interaction effects can never be sufficient to 
cancel the square-root singularity due to the instantaneous elastic response. 



B. The ( — h) and (— ) families 

For the ( — h) family, whose Z function diverges at the slip zone tail only, the analysis 

parallels completely the above one, the $ |_ fields obey a set of equations with exactly the 

same structure as that of Eqs.([29"l-j32|), differing only in the detailed algebraic expressions 
of G(u), C* , D*. The regularization condition analogous to Eq . (|39"D , now to be imposed at 
u — > —a, is again immediately shown to have no solutions. 

Z diverges at both slip edges, hence tow regularization conditions, and one shows 

similarly that the condition obtained from their difference cannot be satisfied. Note, however, 

that the counting argument tells us that, if ( ) patterns could exist, one only of the five 

parameters would be free, i.e. there could exist at most one dynamical pattern at a given 
sliding velocity. 
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C. The (++) family 



As Z ++ vanishes at both slip zone edges, no regularity condition has to be imposed |20 . 
(++) solutions, if any, form a two parameter family.. The analysis of Appendix B again 
leads to an expression of $ ++ with the same structure as Eqs. ( |2"D| - |3"2"|) . One can then write 
explicitly the self consistency equation for Vo- We will skip here the corresponding 
tedious but straightforward algebra, and only quote the final form of Eq. fl2"5D , which can be 
written as : 



Va 



1 + ^ 



t* 1 a, 2 a 

: — (1 — COS — ) COS — 

2/i ckr y 2 J 2 



(41) 



For the systems we are interested in, as already mentioned, /loo/Vo << 1- Then, again, 
under the stick restriction which imposes that r* < 0, condition (|29|) cannot be fulfilled. 



IV. DISCUSSION 

The above analysis leads us to a strong statement, which seems to contradict existing 
qualitative observations. Namely, an interface with Coulomb friction between a viscoelastic 
and a non-deformable material cannot sustain slow sliding via a periodic set of alternating 
non inertial slip pulses and stick regions. 

We believe that the reason for this contradiction must be traced to the fact that the 
Coulomb model of friction which we have assumed misses some physical elements which 
probably play a crucial role in the dynamics of patterns with fracture-like singularities. This 
may appear more clearly when one notices that this model, which describes the interface as 
infinitely rigid below the friction threshold, then, once this is reached, sliding under constant 
stress, is the exact 2D equivalent of the Hill model of bulk ideal plasticity |H| , well known 
to generate numerous artefact instabilities due to its highly singular character. 

Clearly, the main weakness of the Coulomb model lies in its overschematic description of 
the transition between stick and slip i.e., for our patterns, in the details of interface boundary 
conditions at the edges of a slip zone. In the case of the analogous mode-I problem - the 
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Griffith crack - it is well known that discontinuous boundary conditions (vanishing normal 
stresses and displacements in, respectively, the cracked and uncracked regions of the crack 
plane) miss a major physical ingredient, namely the finite range of atomic decohesion and its 
associated energy cost. Taking this into account regularizes the stress field at the fracture 
head, by smearing its square root singularity over the Dugdale-Barenblatt cohesive zone fZ2 



In analogy with this, and based upon the nature of the stress-strain characteristics of 
overconsolidated clays, Palmer and Rice |23f proposed a model for sliding along a concen- 



trated slip surface in which the sliding shear stress is assumed to decrease with relative 
displacement as shown on Figure 2. This enabled them to analyze the "mode-II fracture" 
problem of shear band propagation in such materials. 

Let us assume for the moment that we can modify our Coulomb model in a similar 
manner. That is, let us assume that the shear stress in the sliding state is given by : 

T V 2 = -fT22 + 5T [8(x l ,t)] (42) 

where 8tq is maximum for 8 = and has a small range 8q << 21. We define 8(xi,t) as the 
displacement at the interface point x± from its position when it was in the preceding stick 
zone, i.e. up till the head of the slip zone under consideration reached it. So : 

8(x u t) = u x {x u t) -u 1 {x 1 ,t- = / Q drj V(rf) = 8( V ) (43) 

C J T] 

In Eq . (|21~1) for the velocity field, r* must now be substituted by r* — 8r [8(ri)]. 

In order to fix ideas, let us concentrate on (H — ) solutions. Repeating the analysis of 
Appendix B leads again to expression (0), with : 

H + 4u) -> F + >)-%^/(a) (44) 



vrJ y 1 + v 2 ) a-y )Z+Jy) v ' 



+ y 2 )(a-y)z+-(y) 

The regularity condition then becomes : 
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t: < cos » + 2^^-g(») = m (46) 

The l.h.s. of Eq.(|46|) must, as shown above, be negative. From Eq. (|45|) , 1(a) < 0. So, the 
introduction of a "cost for incipient sliding" via 5tq is sufficient to lift the incompatibility 
which we found to hold for the Coulomb friction model. Clearly, the same formal result 
applies for the (-+) and (-) classes. That is, the localized incremental stress 5t q plays a role 
comparable to that of the cohesive stress in mode-I fracture, namely it smoothes out the 
stress singularity by spreading it over a zone of incipient sliding of small but finite extent. 

Once this formal remark has been made, one should however come back to the possible 
physical interpretation of such a modification of the friction model. A decrease in the 
frictional stress with the slip distance is likely to be associated with a change upon sliding 
of the internal structure of the nanometer-thick adhesive interfacial junction. Moreover, in 
order for the peaked structure of Ti 2 (<5) to reproduce itself at each successive slip zone head, 
the structure of the junction must relax non negligibly on the duration Ar st = (A — 21) /c of 
a stick (typically, Ar st lies in the range of seconds). 

Such a scenario is plausible for junctions composed of long molecules - either because a 
molecular layer of lubricant is present or because the junction is formed by molecular tails 
from the sliding material itself. Then, sliding is likely to give rise to a slip weakening of 
friction associated with molecular elongation and restrengthening by structural relaxation 
during stick. These are precisely the physical ingredients invoked to explain the hysteretic 
frictional dynamics observed in a number of boundary lubrication experiments [24], [p5j1 , |26| 



However, inclusion of slip-weakening of dynamical friction is not the only possible im- 
provement on the Coulomb model susceptible to allow for slip pulses. Indeed, a series of 
recent works by Langer et al |27| [28[] on the viscoplasticity of amorphous solids point very 



convincingly towards the crucial importance of a realistic description - of the rate and state 
type - of the gradual cross-over of the mechanical shear response from mainly elastic to 
mainly dissipative. As already pointed out, solid friction along a continuous interface is 
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nothing but 2D interfacial viscoplasticity, to which the bulk analysis should be transpos- 
able. Hence the need for the elaboration of a phenomenolgy which can bridge realistically 
between static and dynamic solid friction. Work in that direction, based upon experimental 
studies of dynamic interfacial shear response, is presently in progress. 

This discussion naturally leads us to emphasize the need for the development of system- 
atic experimental studies of interfacial slip pulses, and the interest which they present. The 
main questions to be elucidated are : 

(i) the precise conditions for frictional sliding to occur in this mode. This includes sys- 
tematic characterization of the bulk viscoelasticity of systems which do exhibit this behavior, 
and qualification of the relevant range of driving velocities Vq. 

(ii) the Vo-dependence of the apparent friction coefficient r* 2 / | | , and the question of 
pattern selection. Namely, is the slip pattern unique for a given v , or, for example, does it 
depend on the lateral size of the sliding block? In other words, does injection at the back 
free edge of the slider play a crucial role in the selection of the pattern wavelength A, or not? 

Further elucidation of these questions would also be of value to shed further light on the 
still largely open question of the physics of shear interfacial fracture. 
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APPENDIX A: 

We briefly sketch here the derivation of eqs.([7|-^|) for a purely elastic system. Let A, \x be 
its Lame coefficients, related to the Young modulus E and to the Poisson ratio v by : 
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E 



A + n 



E 



2(l + i/) ^ 2 (1 + u) (1 - 2z/) 

The elastic displacement u = (mi,m 2 ) obeys the Lame equation : 

pii = (A + ji) 'V.divu + fj.Au 

with p the mass density. The stresses are given by : 



Tij dui du k / 2 x 

+ ~ h (/3 - 2 J 5 ifc dzuu 



where : 



One then sets 



fx dxk dxi 



P 2 



2 (1 - u) 
1 - 2v 



(Al) 



(A2) 



(A3) 



(A4) 



Ui 



{Xi,X2, t) = U* (X 2 ) + U im (^2) e 



imk(x\—ct) 



(A5) 



with u* the displacement field corresponding to uniform sliding under the homogenenous 
stresses r*. Solving equation (|A2[ ) together with the condition of non separation at the 
interface u 2 | I2 =o = 0, one gets : 

k 2 



m (x 1 , x 2 , t) - u* (x 2 ) = Re Y A r , 



m>l 



S + S_ 



+ e 



-ms-X2 



imk(xi—ct) 



(A6) 



w 2 (xi, x 2 , t) — m 2 = Re A n 



ik 



-ms + X2 _|_ g-ms_x 2 



D imk(xi — ct) 



m>l 



with 



si = A; 2 1 



per 



fc 2 (1-^-1 i?es ± >0 



(A7) 



(A8) 



\ A + 2fi J \ p 

Then, with the help of Eq . (|A3|) , and in the incompressible limit A — > 00, one obtains the 
expression for the interface stresses and the interface sliding velocity : 



T~12 \x 2 =0 



r* 2 + LiRe £ mkA m (- - S -) e imk ^ 

m>l / 



(A9) 



17 



^22 \x 2 =0— T22 + f^Re Y] lUlkA 



m>l 



2 ^ S ~ 

s_ k 



(A10) 



v — cRe imkA Ti 



m>l 



1 - 



imk(xi — ct) 



(All) 



Setting : 



(A12) 



and substituting s_ by s m _ = ka m (Eq.( |A3| )) appropriate to the viscoelastic system directly 
yields expressions (pj)-(|TT|). 



APPENDIX B: 



Following [T^j, the singular integral equation ( p4[) , valid in —a < u < a 

a 

jdy 



a . . 

2/ '" / ° (//) + I" dy<S>(y)k(u,y) = F(u) 



tc J u — y 

—a 



where k is given by Eq.(^) and : 



F(u) 



2V u - T* 
1 + u 2 



is equivalent, for the (H — ) family of solutions, to : 



$ + K *k*$ = K * F 



with 



/(</) 



(Bl) 



(B2) 



(B3) 



2/x 7T ^ Z+_ (y) (tt - y) 
Integrating in the complex y-plane along the contour shown on Figure 3, one finds : 

Z+_ (u) C*u + D* 



(B4) 



K*F 



2^o 1 + u 2 



(B5) 



where C*,D* are given by equations fl34|) , d35|). 
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On the other hand : 



(AT*fc*$) (u) 



2fl 7T 



f dz$(z)J^(y,z) 

J —a 



(B6) 



with 



J^(y,z) = - 1 

7T J 



dy 



z +- (y) (« - y) 



;30 



ds 



_ z-T(s) 

y \+zT(s) 



(B7) 



T (s) = tan — 



(B8) 



Once the order of the y and s-integrals on the r.h.s. of eq. ( |B7| ) has been interchanged, the 
y-integration can be performed explicitly. However, care must be exercised when performing 
this interchange, due to the presence of the two principal values. One uses the following 
identity, which results from the Poincare-Bertrand theorem ]BJ : 



PV 



x" — x' 



PV 



PV 



1 



PV 



1 



v i — x' 
+ it 2 8 (x" - x') 5 (x" - x) 



PV 



(B9) 



One thus obtains: 



J-+ (y, z ) = f ds v' ( s ) i 



dy 



a Z+- (y) (u -y)(y- ± -^ ZL 



l+zT(s) 



+ Y 



(BIO) 



y = - 2 E 



fi' (s p (z,u)) 



p=—oo 



dD/ds \ Sp Z, 



u 



-6{s p (z,u)) 



where the s„ are the zeros of D = u 



z-T{s) 
l+zT(s) 



, i.e. 



(Bll) 



s P (z, u) — —j^ [(f) {z, u) + pn] 



(B12) 



(z, u) = tan 



z — u 



1 + zu 



71 

< 

2 



[z,u) < 



7T 



(B13) 



From this one gets finally : 
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(K * k * $) (it) = — — — — J ^ dz®{{z) // [s (z, u)} [9 (0 (z, u))+P] 

z +- («) r^M?^' ( s ) 6 i^ 2 - ° 2 ) 



[ a dz*(z)Us t [S Z~ (B14) 



where /3 is defined in eq.([33]), and : 



Straightforward integration then results in eqs. ([29| - |32|) . 
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FIGURES 

FIG. 1. Schematic representation of the sliding system 
FIG. 2. Schematic representation of the slip-weakening friction law. 
FIG. 3. Contour for integrals of type K * F. 
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